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We study the transmission properties of a guide consisting oj a dielectric 
rod with rectangular cross section, surrounded by several dielectrics of 
smaller refractive indices. This guide is suitable for integrated optical 
circuitry because of its size, single-mode operation, mechanical stability, 
simplicity, and precise constmction. 

After making some simplifying assumptions, we solve Maxwell's 
equations in closed form and find, that, because of total internal reflection, 
the guide supports two types of hybrid modes which are essentially of the 
tem kind polarized at right angles. Their attenuations are comparable to 
that of a plane wave traveling in the material of which the rod is made. 

If the refractive indexes are chosen properly, the guide can support only 
the fundamental modes of each family with any aspect ratio of the guide 
cross section. By adding thin lossy layers, the guide presents higher loss to 
one of those modes. As an alternative, the guide can be made to support only 
one of the modes if part of the surrounding dielectrics is made a low im- 
pedance medium. 

Finally, we determine the coupling between parallel guiding rods of 
slightly different sizes and dielectrics; at wavelengths around one micron, 
S-dB directional couplers, a few hundred microns long, can be achieved with 
separations of the guides about the same as their widths (a few microns). 

I. INTRODUCTION 

Proposals have been made for dielectric waveguides capable of 
guiding beams in integrated optical circuits very much as waveguides 
and coaxials are used for microwave circuitry. 1-3 Figure 1 shows the 
basic geometries for these waveguides. The guide is a dielectric rod of 
refractive index n immersed in another dielectric of slightly smaller 
refractive index n(l — A) ; both are in contact with a third dielectric 
which may be air (Fig. la) or a dielectric of refractive index n(l — A), 
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Fig. 1 — Dielectric waveguides for integrated optical circuitry. 

(Fig. lb). These geometries are attractive not only because of sim- 
plicity, precision of construction, and mechanical stability, but also 
because by choosing A small enough, single-mode operation can be 
achieved with transverse dimensions of the guide large compared with 
the free space wavelengths, thus relaxing the tolerance requirements. 

Even though in a real guide the cross section of the guiding rod is 
not exactly rectangular and the boundaries between dielectrics are not 
sharply defined, as in Fig. 1, it is worth finding the characteristics of 
the modes in the idealized structure and the requirements to make it a 
single-mode waveguide. 

Furthermore, directional couplers made by bringing two of those 
guides close together, Fig. 2, may become important circuit compo- 
nents. 1 - 2 In this paper we study the transmission through such a 
coupler; the modes in a single guide result as a particular case, when 
the separation between the two guides is so large that the coupling is 
negligible. Through use of a perturbation technique, we also find the 
coupler properties when the two guides are slightly different. 
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Fig. 2 — Directional couplers. 
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The guiding properties of the rectangular cross section guide im- 
mersed in a single dielectric are compared with those derived through 
computer calculations by Goell. 4 Similarly, the coupling properties of 
two guides of square cross section immersed in a single dielectric are 
compared with those of two guides of circular cross section derived by 
Jones and by Bracey and others. 5 ' In both comparisons agreement is 
quite good. 

II. FORMULATION OF THE BOUNDARY VALUE PROBLEM 

For analysis, we redraw in Fig. 3 the cross section of the coupler 
subdivided in many areas. Nine of the areas have refractive indexes 
rii to r? 5 ; we do not specify the refractive indexes in the six shaded 
areas. The reasons for these choices will become obvious. 

A rigorous solution to this boundary value problem requires a com- 
puter; 4 ' 7 nevertheless, it is possible to introduce a drastic simplification 
which enables one to get a closed form solution. This simplification 
arises from observing that, for well-guided modes, the field decays 
exponentially in regions 2, 3, 4, and 5; therefore, most of the power 
travels in regions 1, a small part travels in regions 2, 3, 4, and 5, and 
even less travels in the six shaded areas. Consequently, only a small 
error should be introduced into the calculation of fields in regions 1 
if one does not properly match the fields along the edges of the shaded 
areas. 

The matching made only along the four sides of regions 1 can be 
achieved assuming simple field distribution. Thus the field components 
in regions 1 vaiy sinusoidally in the x and y direction; those in 2 and 4 
vary sinusoidally along x and exponentially along y; and those in 
regions 3 and 5 vary sinusoidally along y and exponentially along .t. 
The propagation constants k xi , k x2 , and k xi along x in media 1, 2, and 
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Fig. 3 — Coupler cross section subdivided for analysis. 
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4 are identical and independent of y. Similarly, the propagation con- 
stants /c„i, fc„a, and k v5 along y in the regions 1, 3, and 5 are also 
identical and independent of x. 

In the appendix we calculate these propagation constants and find, 
as expected, that all the modes are hybrid and that guidance occurs 
because of total internal reflection. Nevertheless, because of another 
approximation which consists of choosing the refractive indexes n 2 , 
n a , n t , and n 5 slightly smaller than n x , total internal reflection occurs 
only when the plane wavelets that make a mode impinge on the inter- 
faces at grazing angles.* Consequently, the largest field components are 
perpendicular to the axis of propagation; the modes are essentially of the 
tem kind and can be grouped in two families, E z vq and E" pil . The main field 
components of the members of the first family are E x and H v , while those 
of the second are E„ and H x . The subindex p and q indicate the number 
of extrema of the electric or magnetic field in the x and y directions, 
respectively. Naturally, E' n and E\ r are the fundamental modes; we 
concentrate on them as we discuss the transmission properties of 
different structures. 

III. GUIDE IMMERSED IN SEVERAL DIELECTRICS 

The guide immersed in several dielectrics (Fig. 4a) is derived from 
Fig. 3 by choosing 

C = oo. (1) 

It supports a discrete number of guided modes which we group in two 
families E X VQ and E PQ plus a continuum of unguided modes. 8,8 

3.1 The W v% Modes 

The main transverse field components of the E" vq modes are E y and H z . 
They are depicted in solid and broken lines, respectively, in Fig. 4a for 
the fundamental mode E\ x . Within the guiding rod each component 
varies sinusoidally both along x and along y. Outside the guide each 
component decays exponentially. Such functional dependence is given 
in equation (38) and depicted in Fig. 4b. We assume n 2 9* n 3 ^ n t 9* n s ; 
consequently the field distributions are not symmetric with respect to 
the planes x = and y «= 0. In Fig. 5a we assume n 2 = n A and n 3 = ti s ; 
the E" pq modes depicted are either symmetric or antisymmetric with 
respect to the same planes. These modes look similar to those in laser 



* This approximation is not very demanding. Even when th ia 50 percent larger 
than na,rus,ru, and ns, the results are valid. 
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Fig. 4 — Guide immersed in different dielectrics: (a) cross section and (b) field 
distribution of the fundamental mode En". 

cavities with rectangular flat mirrors, but our nomenclature is different. 10 
The subindexes p and q indicate the number of extrema each component 
has within the guide. 

Now we describe these modes quantitatively by reproducing the 
propagation constants found for each medium in Section A.l of the 
appendix. Let us call k, the axial propagation constant and k xp and k y , 
the transverse propagation constants along the x and the y directions, 
respectively, in the yth medium (v = 1, 2, • • • 5). Furthermore, let us call 



k, = kn y = — n. 



(2) 



the propagation constant of a plane wave in a medium of refractive 
index n v and free-space wavelength A. 
According to equations (39) through (52) 



k., = (/Cj k x k v ) 



(3) 
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Fig. 5— (a) Field configuration of E pq " modes, (b) Field configuration of E, q * 
modes. 



in which 
and 
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(4) 



K v — K v \ — /C„3 — /C v5 . \D) 

This means that the fields in media 1, 2, and 4 have the same x 



DIELECTRIC WAVEGUIDE 



2077 



dependence and similarly those in media 1, 3, and 5 have identical y 
dependence. These transverse propagation constants are solutions of 
the transcendental equations: 



k x a = pw — tan ' fc x £ 3 — tan ' fc x £ 5 



-i n, 



-1»4 



in which 
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(7) 

(8) 
(9) 



and 



^2,3.4.5 = 



(10) 



(fc? - ^2.3.4.5)* 2(n? -n5.». 4l ,)*' 
In the transcendental equations (6) and (7) , a and b are the trans- 
verse dimensions of the guiding rod, and the tan -1 functions are to be 
taken in the first quadrant. 

What are the physical meanings of £ 3 , Vs > and A 2,3.4,5 ? The amplitude 

5 4 
of each field component in medium 3 (Fig. 4) decreases exponentially 
along x. It decays by 1/e in a distance £ 3 = 1/ I k z3 | . Similarly £ s , t? 2 , 
and 774 measure the "penetration depths" of the field components in 
media 5, 2, and 4, respectively. 

The meaning of A 2 is the following. Consider a symmetric slab derived 
from Fig. 4 by choosing a = °o and n 2 = n 4 . The maximum thickness 
for which the slab supports only the fundamental mode is A 2 . 

Expressions (3) , (8) , and (9) contain k x and k v , which are solutions of 
the transcendental equations (6) and (7) . These cannot be solved exactly 
in closed form. Nevertheless, for well-guided modes, most of the power 
travels within medium 1, implying 
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It is possible then to solve those transcendental equations in closed, 
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though approximate, form. Their solutions are 



*. = 8*(l + ^ ■'■ 



wa 



/■ = 31. f i _l w ^2 + n\A A 



■wn x b 



(12) 



(13) 



For large a and b, the electrical width, k x a, and the electrical height, 
k y b, of the guide are close to p-n- and q/tt, respectively. 

Substituting equations (12) and (13) inequations (3), (8), and (9), 
we obtain explicit expressions for k z , £ 3 , £ B , rj 2 , and 774 : 
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(14) 
(15) 

(16) 



3.2 The E Z VQ Modes 

Except for the fact that the main transverse components are E x and 
H y , the E z vq modes are qualitatively similar to the E" VQ modes (Fig. 5b) ; 
they differ quantitatively. Distinguishing with bold-face type the symbols 
corresponding to E z pa modes, the axial propagation constant and the 
"penetration depth" in media 2, 3, 4, and 5 are, according to equations 
(60), (63), and (64), 



k, = (fc? -k\- kl) 
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in which k x and k„ are solutions of the transcendental equations 

2 2 

k x a = vir — tan -^k x ? 3 — tan — k x £ 5 

k„& = gx — tan -1 k„n 2 — tan -1 k„n 4 . 
The approximate closed form solutions of these equations are 

n\A z +nlA, w 



k x =^(l + 
a 



Trn,a 



and 



ktf = f ( 1 + i^u 



(20) 
(21) 

(22) 

(23) 



Substituting these expressions in equations (17), (18), and (19), we 
derive the explicit results: 
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If 
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3 



«1, 



these results coincide with those in equations (14), (15), and (16), 
indicating that the E z pa and E" pa modes become degenerate. 

3.3 Examples 

The axial propagation constants k, and k, , given in equations (3) 
and (17) and properly normalized, have been plotted in Figs. 6a through 
k as a function of the normalized height of the guide 

6 26 . 2 2\* 



2080 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1969 
I. Hi 




0.4 0.8 



1.2 1.6 2.0 2.4 2.8 3.2 3.6 

b _ 2b / n 2 n 2\ '/2 

A-irr _n4 J 



Fig. 6 — Propagation constant for different modes and guides. 

scendental equation solutions; closed form solutions; 

computer solutions of the boundary value problem. 



tran- 

— Goell's 



for several geometries and surrounding media.* The ordinate in each of 
these figures is 



it varies between and 1. It is when k t = k A , that is, when the guide 

* In these figures we use the same symbol k, for both the E pn " and the E,,,, 1 
modes. 
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is so small that the mode under consideration becomes unguided or, in 
other words, the "penetration depth" in medium 4 is » . It is 1 when the 
guide is so large that k z = k, , which means that all the field travels 
within the guiding rod and the "penetration depths" in media 2, 3, 4, 
and 5 are zero. 

The solid curves have been obtained using the exact numerical 
solutions of the transcendental equations (6), (7), (20), and (21); 
for the transverse propagation constants k x and k„ ; the dashed lines have 
been derived using the closed form approximations (12), (13), (22), 
and (23). In Figs. 6a, 6b, 6e, and 6f, for comparison, we have also 
included the dotted-dashed lines which are the results obtained by 
Goell as computer solutions of the boundary value problem. 4 

The three solutions coincide even for moderately large values of b. 
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Thus, for a guide and mode for which 



k\ - kl 



^ 0.5, 



the closed form approximation is within a few percent of the exact value. 
This gives us confidence to use our results in guides with an aspect 
ratio a/b > 2, in guides surrounded by several dielectrics and in direct- 
ional couplers for which there are no computer calculations available. 
The largest discrepancy between our results and Goell's occurs for 
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and especially for the fundamental modes E x u and E" t . Our approximate 
theory is incapable of predicting the fact that these modes remain 
guided no matter how small the guide's cross section. 

Figures 6a through d cover the cases of rectangular guides totally 
embedded in a single dielectric of slightly lower refractive index. For all 
practical purposes, given p and q, the E' pq and E PQ modes are degenerate, 
and the square cross section provides the widest separation between 
modes. 

Figures 6e through g also consider rectangular guides embedded in a 
single dielectric, but the external refractive index is 1.5 times smaller 
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than the internal one. A glass rod immersed in air is an example. The 
substantial difference of refractive indexes breaks the degeneracy for any 
rectangular cross section. Rectangular waveguides as in Fig. la, with 
three sides in contact with slightly lower refractive indexes and the 
fourth side in contact with air, are covered in Fig. 6h through k. 

The approximate dispersion relation (14) for E U PQ modes, in a rectang- 
ular guide surrounded by four different dielectrics, has been put in 
graphical form in Fig. 7 by plotting the equivalent equation 
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Pig, 7 — Nomograph to dimension a guide immersed in several dielectrics in 
such a way that it supports any prescribed number of modes. 
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The curves plotted for different values of p and q are straight lines 
(solid lines) ; since the values of X and Y are physically meaningful 
when they are positive, the plots are kept within the first quadrant. 

In Fig. 7 the dotted lines depict the equation 



= C. (30) 



Given any guide, we can calculate C which is a function of the dimen- 
sions, refractive indexes, and wavelength. The corresponding dotted 
line intersects all the solid lines representing the different modes. The 
abscissa or ordinate of each intersection yields, after some algebra, the 
propagation constant k~ of each particular mode. If the resulting k~ is 
smaller than the smallest k v , that mode is not guided. 

Another way of using the graph is this: Suppose one wants a guide 
with such dimensions that at a given wavelength only the E^ mode is 
supported. Picking k, = /c„ min , any combination of n t , n 2 , n 3 , n 4 , n B , 
a, and b represented by a point within the triangle limited by the solid 
lines E\ x , E\ 2 , and E\ x will satisfy the proposed single-mode require- 
ment. 

In the graph it is enough to substitute a by 6 and everything we said 
about E V PQ modes is applicable to E x pq modes. 

Figures Ga through k have been used to determine dimensions for 
several guides. All of them have the maximum dimensions compatible 
with exclusive guidance of the E xl and E" t modes. The results are 
collected in Table I. 

In general, the geometry with n 2 < n 4 requires a larger waveguide 
cross section than with n 2 = n 4 . This means reducing the refractive index 
on one side of the guide reduces its ability to guide. The explanation of 
this paradox is found in the known fact that a symmetric slab indeed 
guides "better" than an asymmetric one. Comparing, for example, 
Figs. 6d and 6k, in which the solid curves have been drawn solving 
Maxwell's equations exactly, the E x Pl and E" vl modes can be guided by the 
symmetric slab (Fig. 6d) no matter how small the thickness 6; there is a 
minimum thickness required for the asymmetric slab (Fig. 6k) to guide 
the same modes. 9 

Consider the guide immersed in a single dielectric. In general, the 
guide's height b is inversely proportional to 

1 

{n\-n\f 
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Table I — Typical Dimensions for Several Guides* 
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* Dimensions are for guides capable of supporting only the fundamental modes 
E n x and E n *. 

* All numbers in the table must be multiplied by X/ni. 

For 7i, = 1.5, n 4 = 1, and X = 1m, the largest guide height corresponds 
to the square cross section, and b = a = 0.61/z- This dimension may be 
too small and difficult to control. The tolerance requirements may be 
relaxed by choosing n x — n A « 1 . Nevertheless, this difference cannot be 
made arbitrarily small because the guide loses its ability to negotiate 
sharp bends." 

In all these examples the fundamental modes E z n and E" n are almost 
degenerate, so symmetry imperfections of the guide tend to couple 
these modes. A lossy layer, added to one of the interfaces between 
guiding rod and surrounding dielectrics, should attenuate the mode with 
polarization parallel to that interface. As an alternative, the guide can 
be made to support only the fundamental mode 2?,, by substituting 
medium 2 with a low impedance medium such as a dielectric with large 
refractive index or a metal. 

An example of such a guide and the propagation constant of its modes 
are shown in Fig. 8. By choosing 



a < 



0.7X 

5! - ^)* 



only the E v n mode is guided. If the metal is not perfect, there is power 
leakage into the low impedance medium. The smaller that impedance, 
the smaller the leakage. 

Guides for integrated optics may be easier to build with a/b » 1. We 
can use Fig. 7 to design a guide of arbitrary dimensions a and b which is 
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Fig. 8 — Propagation constant for modes in a guide surrounded by metal and 

dielectrics. transcendental equation solutions; closed 

form solutions; — • — • — Goell's computer solutions of the boundary value 
problem. 

still capable of supporting only the E x u and E\ x modes. An as example, 
let us calculate what the values 

n s = n B = nj(l + A) and n 2 = n 4 = n,(l + A') 
should be, assuming 



Choosing 



A, A' « 1, and r = 5. 



(f 1 - f - «• 



(31) 



one derives from Fig. 7 



C - IgJ :>•;. 



The curve corresponding to C = 25 has been plotted as a dotted line 
in Fig. 7. It intercepts the E\ x line at 

In this expression, by making 

k t = fcn,(l — A), 
the guide supports only the E" u and E z n modes; its height is then 
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b = 1.66- rrrrr- 

n,(A y 

We can choose b arbitrarily by the proper selection of A'- 
For 

X = 1/* n x = 1.5, and b — 5/x, 
from equations (31) and (32) we obtain 

a = 25m, A = 0.002, and A' = 0.05. 



(32) 



IV. DIRECTIONAL COUPLER 

In general, the directional coupler can transmit E x pg and E PQ modes; 
but if the sides a and b of the guides are selected small enough, only the 
fundamental modes E x u and E\ x are guided. Let us concentrate on the 
E v xx mode. The coupler guides two kinds of E\ x modes: one is symmetric 
(Fig. 9c) while the other is antisymmetric (Fig. 9d). Both are essen- 
tially tem modes with main field components E y and H x . The electric 
and magnetic field intensity profiles for both modes are depicted 
qualitatively in Figs. 9b, c, and d. 

Ignoring the small effects introduced by the loose coupling, the 
electrical width k x a and height k y b of each guide, as well as the field 
penetrations £ 3 and r) 2 , coincide with those of the guide described in 

5 4 

.Section III. Similar reasoning applies to the E\ x mode. 

The coupling coefficient K between the two guides and the length L 
necessary for complete transfer of power from one to the other are, 
according to equations (56) and (59), 12 



-iK = 



JL = o fc ' & exp (-c/g B ) 
2L k, a 



(33) 



i + m 

For E V VQ modes, k, and £ 5 are given in equations (3) and (8), and k z is the 
solution of equation (6). Similarly, for E z pq modes, k t , £ 5 , and k x are 
obtained from equations (17), (18), and (20). As expected, the coupling 
decreases exponentially with the ratio c/£ 5 between the guide's separation 
and the field penetration in medium 5. 
The normalized coupling coefficient 
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(34) 
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Fig. 9 — Directional coupler immersed in several dielectrics: (a) cross section, 
(b), (c), and (d) field distributions. 

derived from equation (33) by substituting £ fi for its value given in 
equation (8) has been plotted in Fig. 10 for the E\ q mode, assuming 
n 3 = n 5 and n x ln h is arbitrary. The solid and dotted lines were obtained 
using the exact solution of (6) and the approximate expression (12), 
respectively, for k x . Both sets of curves are close to each other, espec- 
ially for 2a/\(n\ - n*) § ^ 1. 

The dashed-dotted lines are the couplings obtained by A. L. Jones 8 
for two parallel cylinders of refractive index rij = 1.8 embedded in a 
medium n B = 1.5. 5 As expected, if the diameters of the round guides are 
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equal to the widths of the rectangular guides, and if the separations are 
the same, the coupling between the round guides should be slightly 
smaller than that between the rectangular ones. 

The normalized coupling equation (34) for the E\ a mode has been 
plotted in Fig. 11, using for k x the exact solution of equation (20) . 
For rii/wo close to unity, the hues get close to the solid curves in Fig. 10 
as the E" lQ and E\ a modes approach degeneracy. The influence of the 
height b of the guides, the refractive indices n 2 and n 4 , and the value of q 
in the coupling of either mode is not important since they only affect k, . 

To work some examples, assume 



n, = 1.5, 



n 2 = n 3 = n 4 = n 5 = 



1.5 
1.01 ' 



and a = 2b. 




Fig. 10 — Coupling coefficient for Ei q " modes. : wujyuu^ 

from trancendental equations; closed form approximations; 

coupling between two cylindrical rods (A. L. Jones 5 ). 



coupling calculated 
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Fig. 11 — Coupling coefficient for Eu, x modes. 
ni/ns = 1.5; — — — Eif coupling for m/n, = 1.1. 



Eiq" coupling for 



To insure that each guide only supports the i?*, and E{ x modes, the 
normalized dimension b according to Fig. 6b must be chosen to be 



2b 



,2\§ _ a 7K 



y (n? - n 2 4 y = 0.75. 



Consequently 



fc, 



b = 1.77X, a = 3.54X, and jl* = 1- 



From Fig. 10 we obtain the coupler length L for complete power trans- 
fer: 

L = 6540X for c = a and L = 262X for c = -■ 
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How far apart should two guides of length I be spaced to have small 
coupling? If the transfer coefficient | T \ = I \ K | <8C 1, from equation 
(33) we derive 

For the same guide dimensions of the previous example and for 
I = 1 cm, A = 1m, and T = 0.01, 

we derive, from either equation (35) or Fig. 10, that c/a = 2.5. Con- 
sequently, both guides 3.54/x wide and 1 cm long would couple — 40dB 
if their separtion is 8.9/x. 

Now we evaluate how a small change of the refractive index between 
the guides modifies their coupling. Such would be the case if the 
medium between the guides is, for example, an electrooptic material 
and we change the applied field to modulate or switch the output. 

For E z u and E\ x modes, assuming well-guided modes {k x A 6 /ir <<C 1) 
and w, — n 5 /w, <5C 1, the ratio between couplings for two values of 
refractive index in medium 5 (for example, n 5 and n 5 (l + 5)), result from 
equations (34) and (12): 

/~,2 \-i „ s r /o „ Will 

(36) 



K 2 




= exp j- 


-( 


i 


- 


')■' 


That ratio is 


1/2 if 














8 = 0.22 


(n\ 

tfiS 


— 


>) 


t 



>-(M 



(37) 



A directional coupler with coupling coefficient K\ and length L = 
tt/I 2/Ci | would transfer all the power from one guide to the other. 
If the refractive index of the medium between the guides was changed 
from n 5 to n5(l+S) such that equation (37) is satisfied, the power 
would emerge at the end of the input guide. The larger the separation 
c of the guides, and the smaller the difference of refractive indexes 
?'i — ^5, the smaller the change of refractive index required. 

Following the example above, for 

i r 1-5 

n, = 1.5, n 2 = n 3 = n, = n 5 = r-jrr , 

a = 1.5-4 a = 3.54X, and c = a, 
the percentage change of index required is only 8 = 0.0033. 
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V. DIRECTIONAL COUPLER MADE WITH SLIGHTLY DIFFERENT GUIDES 

Consider the directional coupler of Fig. 12 in which the two guides 
have slightly different heights: one measures b + h and the other 
b - h. 

Let us qualitatively plot the coupling coefficient as a function of h, 
Fig. 13. Because of simple arguments of symmetry, the absolute value 
of coupling coefficient is stationary (first derivative zero) around 
h = 0. Therefore, the coupling coefficient between two guides of height 
&i and b 2 is the same as that of the coupling between two identical 
guides of height 1/2 (bi + b 2 ) , provided that | bi - b 2 \ is small 
enough. 

This reasoning applies to guides with different widths, heights, and 
refractive indices, provided that the differences are small enough. Un- 
fortunately, as in most perturbation analysis, we don't know what 
"small enough" is unless we calculate the next higher order term. 

VI. SUMMARY AND CONCLUSIONS 

A dielectric rod (Fig. 4a) of rectangular cross section a by & surrounded 
by different dielectrics supports, through total internal reflection, two 
families of hybrid modes. They are essentially tem modes polarized 
either in the x or the y direction; we call them E x vq and E" vq . The sub- 
indices state the number of extrema (p in the x direction and q in the 
y direction) of the magnetic or electric transverse field components. 

Dispersion curves for guides of different proportions and different 
surrounding dielectric are plotted in Figs. 6a through k. Typical di- 
mensions for several guides capable of supporting only the fundamental 
modes E\ x and E\ x are contained in Table I. 

By picking dielectrics with similar indexes, the guide dimensions can 
be made large compared with X, thus reducing the tolerance require- 
ments. The dimensions a and 6 can be picked arbitrarily and still achieve 




Fig. 12 — Directional coupler with guides of different heights. 
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Ikl 



Fig. 13 — Qualitative behavior of the coupling coefficient as a function of h. 

a guide which supports only the fundamental modes if one can choose 
the refractive indexes. The design is achieved with the help of either 
equation (14) or Fig. 7. 

The penalty one pays with most of these guides is that the funda- 
mental modes are almost degenerate; consequently, symmetry imper- 
fections tend to couple them. A lossy layer added to the interface y = 6/2 
(Fig. 4a) should attenuate the E z n mode more than the E" n . As an 
alternative, the guide can be made to support only the 2UJ, mode by 
metalizing the same interface. Dispersion curves are shown in Fig. 8. 

Since the field is not confined, there is coupling between two of these 
guides (Fig. 3) . Design curves for directional couplers are given in Figs. 
10 and 11. 

Typically, for n x = 1.5, n 2 = n 3 = n 4 = n 5 = 1.5/1.01, a = 3.54A, 
b = a/2 = 1.77 A, and c = a/4 = O.SSA, according to equation (33) the 
length necessary for 3dB coupling is L/2 = 131 A. This length increases 
exponentially with the separation between the guides. 

Increasing the refractive index between the guides by a 3 per 
thousand doubles the coupling. 

What is a reasonable separation to prevent coupling? Using the 
numbers of the previous example, two parallel guides 1 cm long sepa- 
rated by 2.5 times the width of each guide have a coupling of — 40 dB. 

The dielectric waveguides and the directional couplers described 
show great promise as basic elements for integrated optical cir- 
cuitry because they: 

(i) Can be made single mode even though their transverse dimen- 
sions can be large compared with the free space wavelength of opera- 
tion. Consequently, the tolerance requirements can be relaxed. 

(ii) Permit the building of compact optical components. 

(Hi) Are mechanically stable and alignment problems are mini- 
mized. 

(iv) Are relatively simple structures and lend themselves to being 
fabricated with high precision integrated circuit techniques. 

(v) Can include active devices of comparable small dimensions. 
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APPENDIX A 

\ 

Field Analysis of the Directional Coupler 

We solve Maxwell's equations for the directional coupler whose 
cross section is depicted in Fig. 3. The structure is symmetric with 
respect to the x = plane; therefore, the modes have electric fields 
which are either symmetric or antisymmetric with respect to that 
plane. Consequently, the guide we have to study is simpler (Fig 14) : 
if the plane x = is an electric short circuit, the modes of the coupler 
propagating along z are antisymmetric; if the plane x = is a mag- 
netic short circuit, the modes are symmetric. As is known, it is the 
interaction of these symmetric and antisymmetric modes traveling with 
different phase velocities along z that represents the effect of coupling. 

As discussed in Section II, by neglecting the power propagating 
through the shaded areas, the fields must be matched only along the 
sides of region 1. We find that two families of modes can satisfy the 
boundary conditions; we call them E z vq and E V VQ . Each mode in the first 
family has most of its electric field polarized in the x direction, while each 
mode of the second family has the electric field almost completely 
polarized in the y direction. The subindexes p and q characterize the 
member of the family by the number of extrema that these transverse 
field components have along the x and y directions, respectively. For 
example, the E' n mode has its electric field virtually along x, its magnetic 
field along y; the amplitudes of the field have one maximum in each 
direction. 

Each family of modes will be studied separately. 

a.i E v vq Modes: Polarization Along y 

The field components in the ^th of the five areas in Fig. 14 are: 13 

\Mi cos (k x x + a) cos (k u y + /3) for v = 1 

M 2 cos {k x x + a) exp {—ik v2 y) for v = 2 

M 3 cos (k u y + /3) exp (—ik x3 x) for v = 3 

M 4 cos (k x x + a) exp (ik v4 y) for v = 4 

M 6 cos (k v y + j8) sin (k x5 x + 7) for v = 5 



II x , = exp (—ik.z + io)t) 



H vy = 

h„ = -4 d ' H - 



k, dx dy 



(38) 
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Fig. 14 — Coupler cross section with plane x = either an electric or magnetic 
short circuit. 



E„ = 



E„ = 



1 d 2 H xv 
(i)en 2 y k z dx dy 

k 2 nl - kl 



E m . = 



i 



dH T 



H r . 



dy 



in which M, determines the amplitude of the field in the ^th medium; 
a and /3 locate the field maxima and minima in region 1; y equal to 0° 
or 90° implies that the plane x = is an electric (antisymmetric mode) 
or magnetic (symmetric mode) short circuit, respectively; cc is the 
angular frequency; e and n (appearing in k 2 = w 2 e/x) are the permittivity 
and permeability of free space. 

In the vth medium the refractive index is n, , and the propagation 
constants k xv , k yv , and k z are related by 

kl + kl + k 2 = »VS = tf • (39) 

To match the fields at the boundaries between the region 1 and the re- 
gions 2 and 4, we have assumed in equation (38) 

Jc xl = k l2 = k xi = k x (40) 

and similarly to match the fields between media 1, 3, and 5, 

k„i = fc v3 = k u5 = k v . (41) 

Before finding the characteristic equations, let us assume the re- 
fractive index ni of the guide to be slightly larger than the others. 
That is 



'- 1 - 1«1. 



(42) 
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As a consequence only modes made of plane wavelets impinging at 
grazing angles on the surface of medium 1 are guided. Since this im- 
plies that 



K 



«/c, 



(43) 



the field components E x in equation (38) can be neglected. 

Now we match the remaining tangential components along the 
edges of region 1 and from equation (38) we obtain 



7 \ 2 'C» 

tan(/c I/ -±/3j=z4^ 



(44) 



tan 



K 


c 

2 

, c 


+ a 


k x 




a + g 







ictn\ k x5 - + 7 



1 



(45) 



Where there are two choices, the upper ones go together and the lower 
ones go together. 

T. Li pointed out that each of these equations considered separately 
is the characteristic equation of a boundary value problem simpler 
than that of Fig. 14. 8 ' ° Thus for a dielectric slab infinite in the x and z 
directions and refractive indexes as depicted in Fig. 15a, the char- 
acteristic equation for modes with no H y component coincides with 




(a) 




Fig. 15 — Dielectric slabs. 
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equation (44) . Similarly, for two slabs infinite in the y and z directions 
and limited at x = by an electric or magnetic short as in Fig. 15b, 
the characteristic equation for modes with E tt = is equation (45) . 

A similar technique has been used by Schlosser and Unger to find 
the transmission properties of a rectangular dielectric guide immersed 
in another dielectric. 7 If the two guiding rods are so far apart that the 
coupling between them is a perturbation, then 

I k x ,c | » 1 (46) 

and we can rewrite the characteristic equations (44) and (45) with 
the help of equations (39) and (46) , making a and b explicit, as 



_, n; 



k v b = qw — tan ' ^§ k u ij 2 ~ tan 1 -ik v rj 4 ' 

7b\ 7b\ 



lC x d h x ()Cl 



1 I _ 1 I 7.2 (.2 



a 1 + kUi 



where k^ is the solution of 



k x0 a = pir — tan 1 k x0 ^ 3 — tan 1 k xQ ^ s , 
1 1 



V2 = 



n*„ 



-4, 



2 

LL 4 



-kl 



(47) 
(48) 

(49) 
(50) 



and 





r ~\ 
T 

I 5> 


2 1* 

— 1c 2 

"'xQ 



(51) 



^2.3.4.5 - 775 _ .2 7} - / 2 _ 2 U' (52) 

\n>i 1^2,3.4.5) Vn "2.3,4,5/ 

In the transcendental equations (47) to (49) , p and q are the arbi- 
trary integers characterizing the order of the propagating mode, and 
the tan -1 functions are to be taken in the first quadrant. The angles 
ka-a and k y b measure the phase shift of any field component across 
the guiding rod in the x and y directions respectively, or in other words, 
the electrical width and height of each guide of the coupler. On the 
other hand, k x0 a is the electrical width of each guide assuming no inter- 
action between the guides, that is assuming c -» oo. 
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Let us find the physical significance of 772,4 and £ 3i5 . The amplitude 
of each field component in medium 2 (Fig 14) decreases exponentially 
along y. It decays by 1/e in a distance -q-, given by equation (50). 
Similarly 774, £3, and £ 5 measure the "penetration depth" of the field 
components in media 4, 3, and 5, respectively. 

The propagation constant along z for each mode of the coupler is, 
according to equations (39), (40) , and (41) , 

k z = (kl -kl- kl)K (53) 

With the help of equation (48) , the slightly different propagation con- 
stants of the symmetric (y = 90°) and antisymmetric modes (7 = 0) 
are 



M 

*J 

In this expression 



— fczO 



1 1 o Ma Is exp (-c/& 

I ± ^ i2 -1 1 ;„2 y 2 

k z0 a 1 + fcx055 _ 



(54) 



Ko = (k\ - klo - J©* (55) 

is the propagation constant of the E" pq mode of a single guide (c — ■> 00 ) . 
The coupling coefficient K between the two guides and the length L 
necessary for complete transfer of power from one to the other are 
related to the propagation constants k zs and k ia by 



~ iK = fL = 



k„ - k, 



klo £5 exp (— c/g 5 



= 2 ^h^ 



k l0 a 1 + ksotn 



7r ak !0 



1 - 



'ferfA, 



exp 



7TC 



1 - 



/Cjo-^*.^ 



(56) 



As expected, the coupling increases exponentially both by decreasing c 
and by increasing the penetration depth £ 5 in medium 5. 

All these formulas contain either fc^o or k y , which are solutions of 
the transcendental equations (47) and (49). For well-guided modes, 
most of the power travels within medium 1 and consequently 



"'lO-^S 



«1 



(57) 



and 



K U A' 



« 1 



(58) 
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It is possible then to solve those transcendental equations in a closed 
though approximate form by expanding the tan -1 functions in power of 
those small quantities and keeping the first two terms of the expan- 
sions. The explicit solutions of equations (47), (49), (50), (51), (55), 
and (56) are given in Section III. 

a. 2 E X VQ Modes: Polarization in the x Direction 

The field components and propagation constants can be derived from 
those in Section A.l by changing E to H and m to — e, and vicel versa. 
Except for their polarizations, the E x pa and E v pq modes are very similar 
and have comparable propagation constants. Using boldface type to 
distinguish the symbols corresponding to E X VQ modes, from equations 
(56), (55), (47), (49), (50), and (51), we obtain 



.~ _ JL _ o ^£° Is exp (— c/j; 6 ) 
** " 2L Z k !0 al + M* 



where 



k z o — ("a &xo *■„) 

and k l0 and k„ are solutions of the transcendental equations 
k u b = gir — tan -1 k„n 2 — tan -1 k„n 4 



and 



in which 



n - n 

k Ttl a = pv — tan -1 z^^olz — tan ' ^kao£ s 



n\ 



n. 



n 2 = 





7T 

A 2 

4. 


2 



(59) 



(60) 



(61) 



(62) 



(63) 



and 



& = F 



(64) 



As in Section A.l, the transcendental equations (61) and (62) can 
be solved in closed, though approximate, form provided that 





A 3 

l 5> 


2 \ 

-k° 



k x o^4 ; 



« 1 



(65) 
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and 

k„A 2 

«1. (66) 



The explicit results are given in Section III. 

REFERENCES 

1. Miller, S. E., "Integrated Optics: An Introduction," B.S.T.J., this issue, pp. 

2059-2069. 

2. Schineller, E. R., "Summary of the Development of Optical Waveguides and 

Components," Report 1471, Wheeler Laboratories, Inc., April 1967. 

3. Kaplan, R. A., "Optical Waveguide of Macroscopic Dimension in Single- 

Mode Operation," Proc. IEEE, 51, No. 8 (August 1963), pp. 1144-1145. 

4. Goell, J. E., "A Circular-Harmonic Computer Analysis of Rectangular 

Dielectric Waveguides," B.S.TJ., this issue, pp. 2133 7 2160. 

5. Jones, A. L., "Coupling of Optical Fibers and Scattering in Fibers," J. Opt. 

Soc. Amer., 55, No. 3 (March 1965), pp. 261-271. 

6. Bracey, M. F., Cullen, A. L., Gillespie, E. F. F., and Staniforth, J. A., "Sur- 

face-Wave Research in Sheffield," I.R.E. Trans. Antennas and Propagation, 
AP-7, Special Supplement (December 1959), pp. S219-S225. 

7. Schlosscr, W., and Unger, H. G., "Partially Filled Waveguides and Surface 

Waveguides of Rectangular Cross Section," Advances in Microwaves, New 
York: Academic Press, 1966, pp. 319-387. 

8. Collin, R. E., Field Theory of Guided Waves, New York: McGraw-Hill, 1966, 

pp. 470^77. 

9. Nelson, D. F., and McKenna, J., "Electromagnetic Modes of Anisotropic 

Dielectric Waveguides at p-n Junctions," J. Appl. Phys., 38, No. 10 
(September 1967), pp. 4057-4074. 
10 Fox, A. G., and Li, T., "Resonant Modes in a Maser Interferometer," 
B.S.T.J., 40, No. 2 (March 1961), pp. 453-488. 

11. Marcatili, E. A. J., "Bends in Optical Dielectric Guides," B.S.TJ., this issue, 

pp. 2103-2132. 

12. Miller, S. E., "Coupled Wave Theory and Waveguide Applications," B.S.T.J., 

S3, No. 3 (May 1954), pp. 661-719. 

13. Schelkunoff, S. A., Electromagnetic Waves, New York: D. van Nostrand, 

1943, p. 94. 



